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Short-range correlations in 4 He are investigated using many-body wave functions obtained in the 
no-core shell model. The similarity renormalization group (SRG) is used to evolve the Argonne V8’ 
interaction and the density operators. The effects of short-range correlations are reflected in the two- 
body densities in coordinate space as a function of the distance between two nucleons, or alterna¬ 
tively in in momentum space as function of the relative momentum between two nucleons. The SRG 
transformation is performed in two-body approximation. The importance of missing three-body and 
higher-body contributions is investigated by comparing results obtained for different flow parameters 
and by comparing to exact results with the bare Argonne V8’ interaction obtained in the correlated 
Gaussian approach. 
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1. Introduction 

Short-range correlations in nuclei reflect properties of the nucleon-nucleon interaction at short 
distances or high momentum transfers. JLAB experiments [1] measuring ( e, e'pn ) and (e, e'pp) 
knockout of nucleon pairs at high energies found for example a strong dominance of pn- over pp- 
pairs. This has its origin in the strong tensor force that induces corresponding strong tensor corre¬ 
lations in the nuclear many-body wave function. Such short-range correlations arc however difficult 
to describe in many-body approaches. Jastrow-type correlation functions are used in variational or 
Green’s function Monte-Carlo [2,3] and correlated basis function [4] methods. In these approaches 
the short-range correlations are incorporated into the wave functions explicitly. Another popular ap¬ 
proach is to use unitary transformations to soften the interaction as in the unitary correlation operator 
method (UCOM) [5,6] or the similarity renormalization group (SRG) [7,8], The unitary transfor¬ 
mations arc typically not applied to the wave functions but mapped onto the Hamiltonian. The wave 
functions obtained with such a transformed Hamiltonian will not have the strong short-range cor¬ 
relations as obtained with the original interaction. An important point here is that also all other 
observables should be transformed in the same way as the Hamiltonian. This has been studied in 
detail for the two-body system in [9], However, many low-energy observables arc not really sensitive 
to the short-range behavior of the wave function and can be calculated in good approximation with 
the untransformed operators. For other observables that arc sensitive to short-range correlations like 
momentum distributions it is however essential to use the transformed operators. 

In this contribution we study short-range correlations in 4 He with the Argonne V8’ interac¬ 
tion [10] by analyzing two-body densities in coordinate and momentum space. Our calculations use 
many-body wave functions obtained in the no-core shell model (NCSM) [11] with an SRG trans- 
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formed Hamiltonian. By comparing the calculated two-body densities obtained with the bare and the 
SRG transformed density operators with exact results for the hare interaction obtained with corre¬ 
lated Gaussians [12] we can illustrate how the SRG transformations indeed remove the short-range 
correlations from the wave functions and how the short-range correlations can be recovered by using 
transformed density operators. 

2. Method 

The SRG transformed Hamiltonian H a and the SRG transformation matrix U a are obtained by 
solving on a momentum space grid the SRG flow equation in two-body space for the Hamiltonian 
and the transformation matrix: 

= [^’^1 ’ ^ = -OaVa, la = VP? [f,H a \ . (1) 

In case of the Argonne V8’ interaction relative momenta up to k max = 15 I'm ” 1 arc included. The 
standard meta generator //„ given by the commutator of the kinetic energy and the Hamiltonian is 
employed. By restricting the evolution to the two-body space the effective Hamiltonian is obtained in 
two-body approximation where higher order contributions are omitted 

H a = UlHU a = t + vf ] + ... + yf 11 * f + V f ] . (2) 

This transformation has to be performed for other operators as well. Due to the two-body approxima¬ 
tion the SRG transformation is unitary only in the two-body space. For larger systems the result with 
the SRG transformed Hamiltonian will be different from those obtained with the original Hamilto¬ 
nian. The role of missing higher order contributions can be estimated by analyzing the dependence 
of the results on the value of the SRG flow parameter a. If higher-order contributions are small, there 
should be only a weak a-dependence. 

For the calculation of the two-body densities we first solve the many-body problem with the 
effective Hamiltonian using the NCSM 

Ha Wa) = E a |T ( ,>. (3) 

With the obtained wave functions |T„.) the hare and SRG transformed two-body densities (in two- 
body approximation) can then be calculated: 

Pbare - Wal P l^a), Peff - <^ a \ $1 p Ua |Y a > . (4) 


3. Results 

In Fig. 1 the two-body densities in coordinate space have been calculated for 4 He as a function 
of the distance of the nucleons and for the different spin-isospin channels: 

pfrir) = Yj ~ n - r )?/ • (5) 

i<j 

The NCSM calculations with the SRG evolved interactions arc well converged within the /V max = 16 
model space. However with the bare Argonne V8’ interaction the NCSM calculations are not con¬ 
verged and we use here results obtained with the correlated Gaussian method [12,13]. As expected 
the two-body densities calculated with the bare density operators (top panel) show a very strong a- 
dependence. The two-body densities clearly show the correlation hole induced by the repulsive core 
in the Argonne interaction. For the SRG evolved interactions the correlation hole can no longer be 
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Fig. 1. Two-body densities in coordinate space for the different spin-isospin channels in 4 He calculated with 
NCSM wave functions and an SRG evolved Argonne V8’ interaction with flow parameters of a = 0 (bare), 
a = 0.01 fm 4 , a = 0.04 fm 4 and a = 0.20fm 4 . On the top the two-body densities calculated with bare density 
operators, on the bottom with SRG evolved density operators. 
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Fig. 2. The same as Fig. 1 but for the two-body densities in momentum space. 


observed - the short-ranged repulsive core of the interaction has been transformed away. For the two- 
body densities calculated with the transformed density operators (bottom panel) we would expect the 
results to be only weakly dependent on the SRG flow parameter a if the two-body approximation is 
justified. This is indeed the case when we look at the even channels, especially at short distances. The 
calculations with transformed Hamiltonians and transformed density operators show the same corre¬ 
lation hole as obtained with the bare interaction and bare density operators. A certain a-dependence 
can be seen at distances around 1 fm and in the odd-channels. In 4 He only a small number of pairs 
is found in the odd channels and they reflect three-body correlations [12] that have their origin in the 
two-body tensor force. With increasing flow parameter the tensor force in the effective Hamiltonian 
gets weakened and the population in the odd channels is reduced. The a-dependence in the odd chan¬ 
nels therefore directly reflects the role of omitted three-body forces in the transformed Hamiltonian. 
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The two-body densities in momentum space as a function of the relative momentum between the 
nucleons arc calculated with the density operator 

hf T {k) = 2 d (3) - kj ) - k) ft?/ . (6) 

i<j 

Results arc shown in Fig. 2. The two-body momentum distributions obtained with the hare interac¬ 
tion can be divided into a low-momentum region up to the Fermi momentum of about 1.5 fin 1 and 
a mid- to high-momentum region. The first observation is that the low-momentum region is almost 
independent from the SRG flow parameter. The results obtained with bare and SRG transformed den¬ 
sity operators arc very similar. We can also observe that the results for the two even channels arc very 
similar in the low momentum region. This is in contrast to the high-momentum region. First of all 
the two-body densities calculated with the hare density operators arc here strongly dependent on the 
SRG flow parameter. With increasing flow parameter the interaction gets softer and correspondingly 
the high-momentum components become smaller. For the largest flow parameter there arc essentially 
no nucleon pairs with relative momenta larger than Fermi momentum. This is of course different 
when the SRG transformed density operators arc used. In this case the high-momentum components 
observed with the hare interaction arc reproduced for relative momentum larger than about 3 fin 1 . 
In the intermediate momentum region between 1.5 and 3 fin 1 we find a strong a dependence in the 
S - 0, T - 1 channel and a relatively weak a dependence in the S' = 1, T = 0 channel. Furthermore 
the number of pairs in the deuteron-like S = 1 ,T = 0 channel is larger by about a factor of ten com¬ 
pared to the S = 0, T = 1 channel. This can be understood by the dominance of tensor correlations in 
this momentum region. The strong a-dependence in the S = 0, T = 1 channel in the mid-momentum 
region is again related to three-body correlations. For very high momenta above 3 fin 1 three-body 
correlations are less important and the momentum distribution reflects the two-body short-range cor¬ 
relations coming from the repulsive core of the Argonne interaction. 

The method presented in this contribution allows to study short-range correlations for all kind 
of interactions and also for all nuclei that are within the reach of the NCSM. It is also possible to 
calculate two-body densities not only as a function of the distance or the relative momentum of two 
nucleons but also as a function of the pair position or pair momentum. 
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